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A possible mechanism for the horizontal line nodes in triplet superconductor, Sr2Ru04, is 
proposed. We consider the interlayer Coulomb interaction, as well as the on-site Coulomb re- 
pulsion, between electrons in different bands. In the second order perturbation in the interband 
interaction, the effective interaction becomes dependent on cos resulting in the horizontal 
line nodes. 
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Since the discovery of the superconductivity in 
Sr2Ru04 by Maeno et al.,^-' its unique properties have 
been revealed by many experiments. Experiments on 
Knight shift^'' and elastic neutron scattering'^) show that 
the electrons make the triplet pairs with the d-vector 
parallel to z as predicted by Rice and Sigrist*' in anal- 
ogy with superfluid '^He. The superconductivity breaks 
the time-reversal symmetry, as indicated by /xSR exper- 
iment. The energy gap has line nodes as shown by the 
temperature dependence of specific heat,^) and the relax- 
ation rate in nuclear magnetic resonance (NMR).^) The 
existence of the line nodes is confirmed by other experi- 
ments.®"^^) The absence of the angle dependence of the 
thermal conductivity in the magnetic field within the a-b 
plane^"' shows that the line nodes run horizontally on 
the Fermi surface. 

When we assume that the order parameter depends 
only on kx and ky, the phenomenological theory based 
on the group theory in the point group of D4/1 shows that 
existence of line nodes and the breakdown of the time- 
reversal symmetry are in general incompatible with the 
symmetry and that only accidental vertical line nodes are 
possible. Quasi-two-dimensional nature of the Fermi 
surface, however, should be taken into account. The 
Fermi surface in Sr2Ru04 consists of three cylindrical 
sheets, which are open in the kz direction. In this case 
the order parameter should be expanded in the Fourier 
series rather than Taylor series (s-, p-, d- partial waves) 
in kz. In the presence of the interlayer interaction, the 
horizontal line nodes are possible. 

Triplet superconductivity is shown to be caused when 
spin fluctuation is anisotropic and the Fermi surface is 
one-dimensional. Takimoto^^^' showed that triplet 
superconductivity appears in the three-band Hubbard 
model when the on-site Coulomb interaction between 
electrons in different bands is large. Nomura and Ya- 
mada^^-' studied the three-band Hubbard model in the 
perturbation theory up to third order. By solving the lin- 
earized Eliashberg equation numerically, they obtained 
that the triplet-superconductivity with the vertical line- 
node-like structure in the /? band is stabilized. The verti- 
cal line nodes, however, can be wiped out by the mixing 
of the order parameters compatible with the symmetry at 



T < Tc.^^'^^^ Since they consider only two-dimensional 
model, ^■^"^^-' the possibility of the horizontal line node 
has not been studied. 

The superconductivity with horizontal line nodes are 
studied by assuming the interlayer attractive interac- 
tion.^^' ^^"^"^ Zhitomirsky and Rice^^-' has proposed the 
pair hopping model between bands, which they call in- 
terband proximity effect. In that model the active band 
has full gap on the Fermi surface, while the passive band 
has line nodes. They argued that there exists the pair 
hopping term from the pair at r and r+{a, 0, 0) in the ac- 
tive band to the pair at r and r -I- (a/2, a/2, c/2) in the 
passive band, where a and c are the lattice constants. 
The position of line nodes has not been studied exper- 
imentally yet. We have proposed that the existence of 
the horizontal line nodes in the nested part of the Fermi 
surface can be observed by the inelastic neutron experi- 
ment. 

Recently, Kondo^^^ studied the two-dimensional Hub- 
bard model at r = in the second order perturbation 
theory. He obtained that singlet superconductivity is 
stabilized in the wide range of electron filling. 

In this paper we propose a mechanism for the hori- 
zontal line node in Sr2Ru04 by applying the Kondo's 
approach to the multiband extended Hubbard model. 
We show that the effective interaction resulting in the 
line nodes can be derived from the lowest order in the 
interlayer interband Coulomb interaction. 

The interband Coulomb interaction is written as 

^iB- E EE^iB(fc-fc') 

k,k',k" (y,a' l^l' 
^ ^k' ,l,a^^^^,'^^k" +k-k' ,V ,cy'^k" ,V ,a' 1 (1) 

where / is the band index and 

lu u'\ ^r(^ , ^{^x ~ K) <K ' '^'y) 
vaB \K, — k) — V\ \^-a cos — — cos — — 
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In the above V is the on-site interband interaction and 
aV is the interlayer interband interaction. 

The effective interaction in the second order in Vie 
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p{e') is the density of states on the Fermi surface and 
z(e) is the momentum dependence of the order parame- 
ter on the Fermi surface (A(fc) — Aoz{e), Aoz{e) cos 
or Aoz{e) cos ckz for the singlet superconductivity and 
dz{k) = Aoz{e), Aqz{6) cos or Aoz{e) cosckz for the 
triplet superconductivity with d{k) \\ z). By defining 



Eq. (7) becomes symmetric form as 



(9) 



Fig. 1. Effective interaction in the second order in the interband 
interaction Vib- Thick and thin lines are the Green's functions 
for the band I and I' ^ I, respectively. Dotted lines are Vib- 



is shown in Fig. 1. Note that singlet superconductivity 
and triplet superconductivity have the same form in this 
model. We get 

Ves,i{q) = -xi'iqx,qy) {ViB{q) f + u 



= Vq {q^ ,qy) + 2Vi {q^ , qy ) cos — 
+ 2V2iqx,qy)cos cqz, 



(3) 



where q = k — k\ U is the on-site intraband Coulomb 
interaction, 

2 

Vo[qx,qy) ^ -xi'{qx,qy){i + Y^'^^ 2'^ 

(4) 

Vi{qx,qy) = -axi' {qx,qy) COS COS (5) 



V2{qx,qyj^ —Xl'[qx,qy) COS —COS — 



(6) 



and xi'io) is the susceptibility in the V band. We have 
neglected the kz dependence in Ckv , since the Fermi sur- 
face has little warping in Sr2Ru04. Then xi'iQ) is inde- 
pendent of qz and the qz dependence of Vcs^i{q) comes 
from the interlayer interaction. We have also neglected 
the second order perturbation in U and the Coulomb in- 
teraction between electrons at the nearest sites in the 
same plane. The interaction within the same plane 
changes only Vo(fca;,fcj,) and does not affect Vi{kx,ky) 
and V2{kx,ky) in Eq.(3). Since we are interested in the 
mechanism for the horizontal line nodes, we have ne- 
glected these terms. 

We write the two-dimensional intersection of the Fermi 
surface as a function of 9, i.e. k on the Fermi surface 
is written as {kp{e) cosO, kp{e) sinO, kz). As shown by 
Kondo,^^^ the most stable state at T = is given by the 
solution 



de'Vptie, e')p{e')z{e') = Xz{e) {i o, i, 2) 



with the largest eigenvalue A, where 

VFi{e,e') = Vi{kF{e) cos6l-fcF(6'')cos6'', 



(7) 



de'VFi{e,e')w{e')^ Xw{e) (i = o,i,2) (10) 



where 



VF,{e,e') = VF^{e,e')^ p(e)p{e') (11) 

The similar equation is obtained by maximizing the 
average of the effective interaction on the Fermi surface 
with respect to the order parameter 



5{VFi{e,e')) 
Sw{e) 







(12) 



where 



,y JJdede'VFdo,o>{e)wie') 
^^^^^^'^^^^ Jdeiwie))^ ('^^ 

is the average of the effective interaction for the order 
parameter z(e) on the Fermi surface. Eq. (12) is written 



VF^{e, e')w{e')de' = {VF^{e, e'))w{e) (w) 

Comparing Eq.(14) with Eq. (10), we get 

\^-{VF^{e,e')) (15) 

Using the Fourier expansion for z{e) we get the eigen- 
value problem of the non-Hermite matrix^^^ from Eq. (7). 
We can obtain the eigenvalue problem with Hermite ma- 
trix from Eq. (10), 



yMlf)wv=\wu 



where 



2tt 
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dede'VF^{e,9')e 



and 



w{e) 



E 



(16) 



(17) 



(18) 



In Sr2Ru04, the (3 and the 7 band have the cylindrical 
Fermi surfaces with the cross section of 0.457 and 0.667 
of the Brillouin zone, respectively.^'^-' As a first step, 
we use the simple model that two bands have the same 
energy dispersion. We take the simple model that there 
are two bands which have the same energy 

ei{k)^ei,{k) 
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= — 2t (cos akx + cos aky) — At' cos akx cos aky — /U 

(19) 

with t = 1 and t' = 0.32. 

The representation of the D4/1 point group is given 
by four one-dimensional representations, ai, a2, bi, 62 

and one two-dimensional representation e. In each one- 
dimensional representation the eigenvector satisfies the 
relation, 

71" 

w{6) = ri2w{9 + tt) = i]4^w{0 + — ) 

= r]uw{-9)=r]u'w{-^+^) (20) 

where 772, rj^, rjjj and tju' arc the characters of the tt- 
rotation around the 2;-axis, the ^-rotation around the 
.2-axis, the 7r-rotation around the a;-axis, and the tt- 
rotation around the x = y line, respectively for each one- 
dimensional irreducible representations. The irreducible 
representations ai, 02, bi, and 62 behave as extended 
s-wave, g-wave {kxky{k^ — ky)), k^ — ky, and kxky, re- 
spectively. 

We exclude the constant term wi^o, which should be 
included in the Oi symmetry, in the /s^-indcpcndcnt order 
parameter in order to avoid a large intraband on-site in- 
teraction U. Since p{6) is not constant, the effect of U is 
not completely avoided by removing wq. The intraband 
on-site interaction, however, changes only superconduc- 
tivity with the fcz-independent Oi symmetry which is not 
very important in this study as we show below. Thus we 
use the approximation of removing wq . 

For the two-dimensional representation, w{9) should 
satisfy 

w{e) = -w{d tt) = ±w{-9). (21) 

In the present approximation that only the lowest order 
in the order parameter is considered, two-dimensional 
representation corresponds to two degenerate states for 
± in Eq. (21), which behaves as kx and ky, respectively. 
For example, the e symmetry for Vfi{9,9') has the the 
order parameter represented by 

ck ck 
d;,{k) « Aokx cos Aoky cos (22) 

or 

A(fe) w Aokx sin ^,Aoky sin ^ (23) 

In our approximation triplet and singlet superconductiv- 
ity have the same stability. The actual order parameter 
in the two-dimensional representation is expected to be 



dz{k) « Ao {kx + iky) cos ^ 



or 



A(fe) « Ao {kx + iky) sin 



ckz 



(24) 



(25) 



if the higher order terms in the order parameter is taken 
into account for the superconducting condensation en- 
ergy. 

In Figs. 2-4 we plot the maximum eigenvalues for each 
irreducible representation as a function of the electron 




Fig. 2. Maximum eigenvalues of VpoC^,^') as a function of elec- 
tron density Ue for a = 0.1. 




Fig. 3. Maximum eigenvalues of Vfi{6,6') as a function of elec- 
tron density He for a = 0.1. 



density for a = 0.1. It is seen from Figs. 2-4 that 
the superconductivity with e x cos ^ symmetry (and 



ckz 
2 



e X sin ^ symmetry, which degenerates with e x cos 
in the present model) has the max;imum eigenvalue for 
0-8 >n^> 0.5. 

In Fig. 5 we show the phase diagram in the Uf^-a plane. 
The triplet superconductivity with horizontal line nodes 
(e X cos ^) is stabilized in some region of parameters, 
which may be realized in Sr2Ru04. 

We have shown that the effective interaction due to 
the interband interaction can cause the superconductiv- 
ity with horizontal line nodes. The mechanism for the 
superconductivity proposed in the present paper is the 
extension of the idea by Little, ^''^ who proposed the at- 
tractive interaction between electrons via a side-chain 
Coulomb interaction. In our case the different band plays 
the role of the side-chain. 

The degeneracy between e x cos ^ and e x sin ^ 
can be lifted by several effects such as warping of the 
Fermi surface due to interlayer hopping, intra-band 
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Fig. 4. Maximum eigenvalues of Vp2{0,6') as a function of elec- 
tron density Ue for a = 0.1. 



«c 0.1 




Fig. 5. Phase diagram for the extended interband Hubbard model 
with t = 1 and t' = 0.32. 



interactions, higher order perturbation, anisotropy of the 
spin susceptibiUty, etc. 

The order parameter with horizontal hne-node is sta- 
bihzcd by the relatively small interlayer interaction (a w 
0.1). The anisotropic superconductivity due to the 
second order perturbation in the on-site interaction is 
caused by the variation of x{q) from its average in order 
to avoid the large on-site Coulomb repulsion U."^^^ On 
the other hand, the product of the on-site interaction 
and the interlayer interaction can use the full value of 
the susceptibility to get the superconductivity. This is 
why the superconductivity with horizontal line nodes is 
stabilized by a small interlayer interaction. 

We have shown that the triplet superconductivity with 
horizontal line nodes, i.e. e x cos ^ state, is shown to 



be stabilized in a reasonable parameter range. We have 
neglected the fact that x(,q) is large due to the nesting 
nature of the a and /3 bands. ^^'^^^ If we consider the 
enhancement of xiQ)i the effective interaction in the 7 
band becomes large due to the mechanism discussed in 
this paper. 

The interband proximity effect"'^^^ is compatible with 
the effective interaction studied in this paper. When 
both interband interlayer interaction and interband prox- 
imity effect are considered, the superconductivity with 
the horizontal line nodes will be much favored. 
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